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SUMMARY 


An experiment is suggested vherein measurements can be made on parameters 
leading to the determination of the homogeneous reaction and diffusion rates of 
oxygen atoms . This is done by studying the decay in the periodic time variation 
in concentration of atoms . The system to be investigated includes a ternary gas 
mixture of diatomic molecules, dissociated atoms, and an inert gas. A theoretical 
investigation is made, starting with equations obtained from rigorous kinetic the- 
ory and chemical kinetics, of the unsteady diffusion and recombination processes 
applicable to the experiment. Boundary conditions are investigated. A qualita- 
tive result is given which illustrates the dependence of the decay in the 
concentration variation on sink strength and frequency. 


INTRODUCTION 


The present theoretical investigation considers methods for analyzing data 
from a proposed experiment on the homogeneous recombination of oxygen atoms to 
determine the diffusion and recombination rates of the atoms at room temperature. 
In the proposed experiment a periodic concentration of atoms is imposed by dis- 
sociation at one end of a test cell. The atoms diffuse Into the test cell and 
recombine, thus providing a continual time -varying supply of freshly recombined 
molecules. The variation with time is a device for distinguishing between recom- 
bined molecules and those not immediately Involved in the recombination process. 
However, to interpret the data from such an experiment, a solution of the unsteady 
multicomponent diffusion equations is required. The appropriate equations are 
derived herein and a simplified qualitative solution is given. Although continual 
reference will be made to the experiment motivating this analysis, the equations 
derived in this report, except the time -dependent boundary conditions, are general 
and apply to any ternary diffusion phenomena as well as being unrestricted by the 
geometry of any experiment . 



The problem to be considered is that of unsteady multicomponent diffusion in 
the presence of chemical reaction inside a tube of circular cross section with 
the time -dependent condition of the gases specified at the entrance of the tube • 
The governing diffusion equations for this problem are derived herein starting 
from a set of fundamental relations given by rigorous kinetic theory. Assumptions 
appropriate for the experimental conditions are used to simplify these relations. 
To the lowest order the resulting equation for atom concentration is linear except 
for a sink term, which represents the atom recombination process . Because of the 
nonlinearity introduced by the sink term, a numerical procedure utilizing elec- 
tronic machine computing appears best adapted to the solution of this equation. 

One can linearize the sink term by use of a rough approximation and, from this 
linearized form of the equation, obtain a qualitative understanding of the more 
difficult problem. To illustrate the effects of recombination and diffusion in 
an unsteady one -dimensional problem, an expression for one component of the 
ternary mixture (i.e., the atom concentration) is derived and discussed. This 
expression shows the dependence of the decay of the labile species variation on 
a linearized sink term and on the frequency of the periodic variation. 
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dimensionless binary diffusion constant ratio. 
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dimensionless constant, 1 
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dimensionless constant, 1 
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species gradient relation arising from rigorous kinetic theory (ref. 5, 
P. ^9) 

dissociation time function 

Maxwellian distribution function for species i 
constant appearing in the defining relation for v 
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dimensionless source strength, — — Ki 
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wall reaction rate 
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coordinate vector 
stream velocity 

mean molecular speed for species 
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characteristic diffusion constant, 


aJ 12D13 

f>23 


multicomponent diffusion coefficient 
binary diffusion coefficient 
source or sink term 

characteristic length appearing in dimensionless coordinate 
number of species comprising the multicomponent diffusion system 
a single direction component of "effusive" flux vector for species i 
peculiar velocity, v - vj_ 


reduced velocity. 



dimensionless angular frequency, 


D 


OJ 


wall reaction fraction 

Laplacian coefficient appearing in partial differential equation for v 

sink strength coefficient appearing in partial differential equation 
for v 

dimensionless perturbation parameter (see eq. (10)) 

dimensionless number density, — 

o n 

equivalent dimensionless concentration variable, ? + g °u 3 

homogeneous part of equivalent dimensionless concentration variable 
which appears in solution for v 

particular part of equivalent dimensionless concentration variable which 
appears in solution for v 
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I" dimensionless coordinate variable, £ 

3 

p mass density, ^Tnjmj 

J=i 

t dimensionless time variable, ^ t 

J_i 

w angular frequency 

A dimensionless velocity, ^ v 

ft square of ratio of diffusion velocity to isothermal sound velocity, 

opD 2 
0 nl 2 kT 

a / . * a 

~ spatial gradient (in Cartesian coordinates x — 

Subscripts 

1 species kind 

x atoms 

2 molecules 

3 inert species 

Subscript appearing before a symbol indicates order in terms of an expansion 
vith respect to A. Superscript appearing before a symbol indicates order in 
terms of an expansion with respect to ft. 

ANALYSIS 
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Before proceeding with the detailed analysis, it may be helpful to consider 
the physical concept of the problem. Figure 1 shows schematically the pertinent 
features of the proposed experimental apparatus. (A more detaxled description of 
a similar apparatus may be obtained from ref. 1.) The essential features of the 
experimental technique are as follows: A mixture of oxygen molecules and some 

inert gas is allowed to leak into a streamtube; flow through the tube is main- 
tained by a vacuum pumping system. The molecules in the gas stream are then 
partially dissociated by means of a modulated radio -frequency signal. This pro- 
cess supplies a periodically varying concentration of atoms at the base of a 
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cylindrical diffusion tube wherein the atoms diffuse and recombine. The variation 
in concentration of either atoms or molecules, depending on the detector used, is 
observed. The following analysis will show that it is then possible to determine 
the n sink strength” associated with the recombination process. The advantage of 
studying (for example, by means of a detector utilizing oxygen molecule ultra- 
violet absorption which is described in ref. 2) the variation in concentration is 
that it is possible to distinguish between the species directly involved in the 
recombination process and those not immediately involved. 


Partial Differential Equations Governing Diffusion 


The analysis that follows starts with a set of fundamental equations 
describing diffusion processes in general. Consideration is then given to a 
physically consistent linearization of the derivative terms contained in these 
equations. The nonderivative or sink term is then analyzed in detail and an 
approximate linearization of this term presented. Finally, boundary conditions 
appropriate to a particular problem are introduced and discussed. 

Development of the derivative terms. - The following equations for diffusion 
processes, resulting from rigorous kinetic theory, are discussed by Hirschfelder, 
Curtiss, and Bird in reference 3* With the neglect of body forces, viscous 
effects, and thermal gradients, the fundamental equations are: 

(l) Equation of continuity of individual species 


hi + A 

bt br 



+ n 


±vi) 
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Ki 
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where 
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n ,i m ,h d Zn p 
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( 2 ) 


and where Kq is the total increase of particles of species i per unit volume 
per unit time due to chemical reaction. 


(2) Equation of motion or the momentum equation 


5v + ^.Sv-_l5p 
bt br P br 


( 3 ) The equation of state 


p = nkT 


( 3 ) 


( 4 ) 
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(4) The defining relations for number and mass densities 


N 



i=i 


( 5 ) 


( 6 ) 


( 5 ) The expressions giving the multicomponent diffusion constants in terms of the 
binary diffusion constants for a ternary mixture with molecules i, j , and k 


Du = 0 


E>ij = D ij 


D ij ~ D ji 


1 + 


n k (ffij D ik " D ij) 


n i D jk + n j^ik + n k®ij 


(7a) 

(7b) 

(7c) 


(Eq. ( 7 b) follows from ref. 3, eq. (11.2-4), p. 716 . ) 

(6) Conservation of mass in chemical reaction 



(8) 


Equations (l) through ( 8 ) form a system of coupled nonlinear partial -differential 
equations . When supplemented with the chemical equations governing possible 
chemical reactions, these relations are sufficient to define the three -component 
diffusion problem under consideration. 

The diffusion process is investigated under conditions where the 
concentrations of the various species are forced to vary in a prescribed periodic 
manner at one of the boundaries. The reaction of the species of primary interest 
is assumed to be governed by the chemical equation 


[A2] -* [Ax] + [Ax] 


(9) 


where the subscript 2 labels the undissociated species in the three -component 
system and the subscript 1, the dissociated species. This subscript notation 
will be kept consistent throughout the report (in addition a subscript 3 will be 
introduced which will be associated with a chemically inert species). Equation ( 9 ) 
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indicates that a unit concentration of species 2 dissociates into two unit 
concentrations of species 1. Depending on the gases used, more than one chemical 
equation may be necessary to describe the chemical reactions possible with such 
species as 1 and 2. The analysis will include the more significant of those . 
reactions possible with diatomic molecules, their dissociated atoms, and an inert 
buffer gas. 

Expressions for the quantities Kq appearing in equations (l) and (8) are 
also needed. These are determined by the form of equation ( 9 ) and a consideration 
of chemical kinetics . We shall return to this point later . 

One can not make immediate assumptions as to the constancy of pressure, 
stream velocity, or mass density without quickly running into inconsistencies . 

For this reason only those minimal assumptions necessary are made within the 
framework of a plausable model. It will be assumed that the molecules can be 
dissociated in a periodic manner without external pressure disturbances and that 
the degree of dissociation will be small. Thus the symbol A is assigned to 
designate the ratio 


A = n x(o)max ( 10 ) 

n(o) 

where n-, ( o)-may is the maximum concentration of species 1 and n(o) is the total 
concentration ( see eq. ( 5 )), both evaluated on one boundary of the system (e.g., 
the input end of a cylindrical diffusion tube). It is assumed that A « 1. All 
of the dependent variables are then expanded in terms of A as follows: 


nq = ( G ni) + (ini)A + . . . i = 1,2,3 (lla) 

v = ( G v) + (iv)A + . . . (lib) 

P = (op) + (ip)^ + • • • (lie) 

n = (o n ) + (in)A + • • . (ill) 

P = (oP) + (iP) A + • ■ • ( lle ) 

Kq = ( 0 Kq) + (iKq)A + . . . i = 1,2 (Hf) 

Dqj = (oDij) + (iE>ij)A + . . . i,j = 1,2,3 i f J (Hg) 


Notice that the subscript denoting the power of A is written on the left of each 
expansion coefficient to avoid confusion with the subscripts on the rights which 
refer to the species involved. 

At A = 0, the gas would be in a uniform state. Therefore, all the zero 
order terms are constant; that is, gradients exist only for the higher order terms. 
When A = 0 (no dissociation), species 1 will not be present. It follows that 

0 ni =0 (12) 
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In the experiment it -would be convenient to include a stream velocity with the 
gas moving axially down a cylindrical tube. It is expected that viscous effects 
would lead to a parabolic velocity profile at the low Reynolds numbers involved. 
However, this complication will be ignored for the present and the equations for 
a constant value of 0 v derived. For values of Q v different from zero, the 
results strictly apply only for a flow far from any solid boundaries. 

The expanded variables (eqs. (lla) through (llg)) are now substituted into 
the original basic equations (eqs. (l) through (8)) and coefficients of like 
powers of A are equated. Since the zero-order terms of the variables are con- 
stant, the zero -order terms occurring in equations (l), (2), and (3) vanish. 

When the relevant substitutions are made, the basic equations become: 

(l) Equations of continuity of individual species 


d(ini) + _3_ 

St Sr 


( i n i) 


oV + 


3 



J=i 


oDij 


_s_ 

_Sr 




iK x 

(13) 


d(in 2 ) + _S_ 
Sr 


( 0 n 2 ) jv + (ins) Q v + ( 0 n 2 ) iV 2 


xK a 


(lb) 


K jjg.). + A 
St Sr 


( O n 3 ) l v + ( l n 3 ) O v + (o n 3) 1 V 3 


= 0 


(15) 


vhere 


lVi=2^3 = 


on 2 

(oP) oni 


0 Dij 

<3=i 


_s_ 

dr 




(16) 


(2) Momentum equation 


d(iy) + , -< S( x v) _ kT S( x n) 
St ° Sr oP Sr 


(IT) 


(3) Equation of state 


0 p = o nkT 


x p = inkT 
etc . 
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(18a) 

(l8b) 



( 4 ) Defining relations for number and mass densities 

G n = Q n 2 + o n 3 


(19a) 


in = 



etc . 


(19b) 


0 P = 0 n 2 m 2 + O n 3 m 3 



j 


etc . 

(5) Multicomponent diffusion constants 


(20a) 

(20b) 


0^12 


D 


12 


1 + 


/ m.3 

° n3 Dl3 ” Dl2 
O n 2D 13 + O n 3 Di 2 


a D 12 D 13 o P 

1^2 D 23 O n 


(21a) 


3 El3 = D 13 


1 + 




o n 2 Dl2 “ Dl3 

0112^x3 + 0^3^12 


a p i 2 p i 3 pP 
m3 D23 cP- 


(21b) 


0^21 “ ^21 


1 + 


^3 


O n 3 l JJ7 ^23 - ^12 

\ x 

on2Di3 + o n 3 D 12 


a D2i ( „ r, , „ "to 


on 5^ I 0 " 2013 + ° n3 ST D23 


(21c) 


)I 9 31 - D31 


1 + 


m 2 


O n 2 1 ^32 - D 13 

O n 2^13 + o n 3 ^ 12 


Jr dJ (° n3Dl2 + ° ns 5 T D32 


(21d) 
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0 E 2 3 = ^23 


( 21e) 

oD 32 = D 32 


(21f) 

Dij = “ 1,2*3 i £ 0 


(21g) 

where the dimensionless constant a is defined as 



D 23 0 n 

a = — — 

O n 2^X3 + o^^x^ 


(22) 

(b) Conservation of mass in chemical reaction 



o 

ti 

o 

ro 

11 

o 


(23a) 

m i(iKi) + m 2 (iK^) = 0 


(23b) 

etc . 



It is convenient at this point to make equations (13) through (17) 
dimensionless. For this purpose a characteristic length (L) and a characteristic 
diffusion constant (D) are introduced. In terms of these two constants and the 
zero-order expressions for number density and mass density, all the variables can 
be made dimensionless. The transformation equations for the dimensionless 
variables In terms of the independent variables are: 



(24a) 

"W* 


(24b) 

II 

3 


(24c) 

where W is a dimensionless frequency, which is needed later 
variables become 

. The dependent 

|i = n/ Q n 


(25a) 

V-i = n i/ 0 n i = 1,2,3 


(25b) 

R = p/ 0 p 


(25c) 

A = (L/D)v 


(25a) 

k ± = (L 2 / 0 nD)K^ 


(25e) 
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The subscripts that will be used before and after the symbols representing the 
dimensionless equivalents of the physical variables will carry precisely the same 
significance as they did for the physical variables . 

A parameter ft will now be defined that represents the ratio of the squares 
of the diffusion velocity to isothermal sound velocity; that is, 


a = Qp _D t- 

o n L 2 kT 



( 26 ) 


The significance of this parameter is that if ft is small, then a pressure pulse 
(or acoustic wave) travels at a rate large compared to any speed associated with 
the diffusion process. The limit ft = 0 implies that the speed of sound is 
infinite relative to any speed associated with diffusion. It is interesting to 
note that in the opposite limit of free molecule flow, ft is of order 1. In this 
limit, the fundamental equations do not apply. In that case the velocity asso- 
ciated with diffusion is of the same order of magnitude as the speed of sound. 
Generally speaking, under those conditions where the fundamental equations 
(eqs . (l) through (8)) apply, ft is small compared to 1. 

Transforming the equations of continuity of individual species into 
dimensionless form, equations ( 13 ) through ( 15 ) become, 


(lhi) + o A ’ ^ (lMa) + ^ m 2 ° 7-^ 


ID] 


8t 




oP D 


( 1^2 ) “ m 2 ( 0^2 ) ' Js ( lt-0 

oP 6^ d| 


A _ oE 


o 




• ^=r (aP-3) " Tp ^3(0^2) ' Jjr ( l^ 1 . 

,d| ° p o| d| 


- i k i 

(27) 


+ oA • i 


+ 





23 


D 



6^ 


( 1^3) 


^ 103 ( 0 ^ 3 ) 


dt 



= 1^2 
( 28 ) 
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(1U3) + o A 


( lM-3 ) + (0^3) JL • (lA) 
O | o| 



0D31 3 
D 


j 5 _ 


(in) 


+ 



0D32 

D 




^2 ( oU-2 ) 
o M 


a 

at 


_a_ 

at 



= 0 


(29) 


In dimensionless form, the momentum equation (17) becomes 


a 


— (iA) + (oA) * (iA) 

St bi 


+ (ip) = 0 


( 30 ) 


Since ft is small, equation (30) can be greatly simplified by an expansion in 
terms of 0 as follows: 


iP = (ip) + (ip)ft + - • • ( 3 l) 

etc . 

The upper left superscript is used here to denote the power of ft Involved. 

When these expressions are substituted Into equation (3O), and coefficients of 
each power of ft are equated to zero, the following relations are obtained: 

4 (?n) = 0 ( 3 2a) 

at 

1 ; (?A) + (qA) • JL (°A) + It (in ) = 0 (32b) 

Thus to lowest order in ft, (Jp) is independent of coordinate. If the first- 
order expression for the equation of state (eq. (l8b)) is transformed into dimen- 
sionless form, equation (32a) is equivalent to stating that the pressure is con- 
stant to lowest order in ft. The lowest order approximation indicates that if a 
pressure pulse Is induced in the system, this pulse will carry instantaneously 
throughout the system. If such pulses are not induced, pressure is independent 
of both space and time. From equations (19b) and (25b), the following identity 
results : 


im(fh) + iH 2 (T,t) + lMUlh) = in = constant (33) 


Since the zero -order pressure term in ft is constant and since the expansions 
are about a uniform state, the constant Is set equal to zero. 
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The continuity equation for species 1 (eq. (27) )> when the terms of order 
ft or greater are neglected, becomes 


(1H1) + o A ’ Vf (1^1) 
ot o| 


o£ 

oP 


m 2 


0^12 9 
d af 



+ 11^3 


0 Pl 3 8 

D df 


8 

8| 



= xki 


( 34 ) 


Substituting the multicomponent diffusion constants (eqs. ( 21 a) and ( 21 b)) and 
simplifying, we obtain 



->A 


81 


(&i) 


3 ^ 120 13 8 
DD 23 < 3 1 


8 / o ^ o. 

■ (iMa) - i k i 

8| 


( 35 ) 


It is interesting to observe that equation ( 35 ) is independent of the first-order 
concentrations of other species as a consequence of the assumption that A and 0 , 
are small. Then this equation can be solved for xhi- The continuity equations 
for species 2 and 3 j equations (28) and (29 ) } respectively, become 


A (1H2) + c>A * A (1^2) + (0^2) • 

ot o| 


(?A) 


+ o£ m 0C21 4 
0 p d 8^ 


r= (1M-1) + ^ m3 

8| oP 


iD 


23 


d 


A (1P3) - lk 2 

a^ 

(36) 


A (1^3) + o-^ ' A (1M3) + (0^3) A ' ( 1-h) 

aT ag ai 


, o n _ 0D31 a 
+ mi ~ 5 f 




0D32 a 

at 


D 


• | <a*> - 0 

( 37 ) 


There are four equations (eqs. (33), (35), (36), and ( 37 )) Tor the dimensionless 
dependent variables representing concentration, iPi, i|J.2> 1^3* and stream veloc- 
ity iA. It is convenient to combine equations (36) and ( 37 ) to eliminate 3A 
and thereby reduce the number of equations and number of dependent variables to 
three. To this end equation ( 37 ) is multiplied by the ratio o^s/oM-s and the 
resulting expression subtracted from equation (36). Using equations ( 21 c) 
through (21g), we express the multicomponent diffusion coefficients in terms of 
the binary diffusion coefficients, and substituting equation (33) for (1P3) we 
obtain an equation in (?p 2 ) and (?p 3 ) alone* The resulting equation can be 
expressed as a sum of two terms, one dependent on (1P2) and the other on (1P.3), 
to obtain 
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( lh2 ) 

dT 


-A 




0^12 5 

D idt 



OP-2 

0^3 



+ 0^ 



aDl3 5 
D LSf 



= 2*2 
(38) 


The separated portions of equation (38) are similar to each other and, in 
addition, each is equivalent in form to the left side of equation (35) • This 
suggests that a differential equation in a new variable v can be found that 
will be equivalent in form to equation (35)- The new variable v will be a 
linear combination of the variables (jV 2 ) an & (1^3) • Substitution of equa- 
tion (33) into equation (35) yields 


A (S*) 


+ 0 ^ * 



3^12^13 

D-D 23 


— . JxL 

M ' 9f 



+ 




a-DigDia d 
I®23 3 I 



-iki 

(39) 


To facilitate the following manipulations let the coefficients of the second, 
third, and fourth bracketed expressions in equation (39) be Pi, Qi, and Si, 
respectively, and the coefficients of the identical bracketed terms in equa- 
tions (38) be P2 , Q2 ? and that is, let 


Pi 


g _ ah 12^13 

1 d 23 d 


(40a) 


Qi = 1 


(40b) 


P 2 


aDi2 

D 


(40c) 


Q2 - 


ph-2 
oh 3 


(40d) 


C _ oh2 al>13 

Ss; ■ yj — 


( 40e) 
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An equivalent form for the differential equation in v can be achieved by 
multiplying equation ( 38 ) by an arbitrary constant F and then adding the 
resulting expression to equation ( 39 ) t o obtain 


( 1 ^ 2 ) + cA * (ite) 

OT 


P 2 + FPi 
F + 1 


.1 ■ l 


Q 2 + FQi 
F + 1 


( 1^3 ) + cA * ( 1^3 ) 

OT o£ 


S g + FS; 
F + 1 


• 4 (?n 3 ) 


= F(gk g ) - (iki) 
F + 1 

(4l) 


It then follows that if F is required to satisfy the relation 


Q 2 + FQi _ S 2 + FSi 
F + 1 P 2 + FPi 


(4la) 


the new variable will be defined by 


= (?p 2 ) + fa -. Jfa (o ) 

V F + 1 / 


(42) 


and v will satisfy the differential equation 

Sv -T . 3v _ (Ps + FPA _B_ _ _6_ = (gkg) - F(gki) (4 } 

Bt ° ’ at \ F + 1 ) St ' Sf F + 1 

which on the left side is of the same form as the left side of equation ( 35 ); 

F is determined in terms of the coefficient expressions by solving equation (4la) 
and is given by 


F 


“ P 2 ^2 

P 2 + Pi(Qe - 1) - S 2 


(44) 


The relevant coefficients in equations (42) and (43) and F can be evaluated 
using equations (40a) through (40e) to obtain 


F 


Q 2 + FQj_ 
F + 1 


OfeC 

(45) 

d - 0 P 2 C 


— = g 

(46) 

0P3 \ d / 
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P 2 + FPj. 
F + 1 


aPi2 

D 


(l - o^c) = 


5i 


Here nev constants defined by 


c = 1 


£13 

^12 


(47) 


(48a) 


d = 1 


P X 3 

D 2 3 


(48b) 


are introduced. Using the conservation of mass relation, equation ( 23 b), °kx 
can be expressed in terms of ikg • The problem then reduces to solving the two 
partial differential equations 


(1U1) + cA ' ^ (ldl) 

OT 0 | 


aU 12^13 d 
DD 23 


tpr ( lUx) - l^i 


(49a) 


_a_ 

St 


v + 



aPi 2 

D 


(1 






nil + OM- 2 C j 

(l - sf) 

d 



(49b) 


where 


V = (xh 2 ) + s( xM-3 ) (49c) 

along with appropriate expressions for the source terms and the boundary 
conditions . 

The solutions of equations (49) are achieved in a straightforward manner. 

If the "sink" or "source" term has the following functional dependence 


xkx = xki(2px) 


(50) 


which would be implied in any chemical equation ^ then the problem requires 
solving one partial differential equation (linear or nonlinear depending on the 
functional dependence included in equation ( 50 )) for iPi- Let the solution be 
expressed as 


il^i = iMa(!,T; ci,c 2 ) 


(51) 


l 6 



where cx and c^ are parameters and are the constant coefficients of the 
Laplacian and sink term, that is. 


cx 


&D 12^13 
hD 2 3 


(5ia) 


C2 “ 1 


(51b) 


On substitution of this solution into equation (50) the final problem will then 
be equivalent to solving a linear inhomogeneous differential equation (eq. (49b)) 
with the appropriate boundary conditions. This solution will be a sum of two 
separate parts, that is, 


v(|,t) = Vp ( I , T ) + v h (|,t) (52) 

where the particular solution Vp(|,T) and the homogeneous solution vjj(£,t) are 
solutions of 



<5t 


+ 



D 


(1 




~xki(ihi) 


mi 

m 2 


+ 




£ 

d 




(53a) 


and 


ha 

St 


cA 


sr 


V H 


aD 
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(1 - oM-ac) ^ V H = 0 

D o| o£ 


(53b) 


respectively, with proper boundary conditions. It turns out that a particular 
solution is given by v p = -iPi(T, t; &i,& 2 ) where 5 X and B 2 are different from 
ci and c 2 . The equations to be solved are as follows: 


£ + ^ | - 1 ■ | ^ w 


h (v p ) + 0 A • h (Vp) - &i 

a T ^ ai 






A (v H ) + <4 • | (v H ) 


Bj. A 
85 




= 0 


(54a) 


(54b) 


(54c) 
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where 


D 


aDigDxa 

P 23 


(54d) 


8 i = — (l - ofec) 


1 ~ Q^gC 
1 - d 


(5^e) 


S 2 


Ei + 

ITlg 




(54f) 


The arbitrary characteristic diffusion constant D is chosen according to 
equation (54d) in order to simplify equation (54a) . 

Development of the sink term .- The preceding analysis has been directed 
toward reducing equations (l) through (8) into the more manageable equations (54) 
for solving the multicomponent diffusion problem, where the concentrations of the 
individual components are fluctuating in some prescribed periodic manner at one 
boundary. The effects of the chemical reactions have not yet been included in 
the problem beyond carrying a lumped term Kq which represents the rate of 
change in concentration of species i due to the chemical processes possible 
with the reactants species 1 and 2. The nature of this "sink 11 or "source" term 
and its expanded version Jk^ will now be discussed. 

The first chemical reaction to be considered is recombination, which is 
described by the chemical equation 


[Ai] + [AjJ - 1 - [M] -* [A 2 *] + [M] + excess energy as heat (55) 


The symbol M is used to denote a third body, which may be any one of the gas 
particles in the system, and which is catalytic in its effect on the reaction. 

The asterisk denotes that the diatomic molecule may be in an excited state (elec- 
tronic or vibrational) and the brackets denote "concentration of" as in equa- 
tion ( 9 ) • The excited molecule subsequently decays to its ground state either by 
radiating or through collision as illustrated by the following equations. 

[A 2 *] [A 2 ] + radiated energy ( 56 a) 

[A 2 *] + [M] [A 2 ] + [M] + heat energy ( 56 b) 
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Generally speaking, because of radiation or binary collisions, the lifetime of an 
excited molecule is small compared to any time associated with diffusion. There- 
fore it is reasonable to expect that after recombining, the excited molecule will 
decay to its ground state before diffusing any appreciable distance in the system 
Equation (55) can then be replaced without appreciable error by an equivalent 
equation 


[Ax] + [A x ] + [M] -+ [A 2 ] + [M] + heat of recombination (57) 


The applicability of this simplification is opportune since it eliminates the 
necessity for solving a four component diffusion problem. There may be apparent 
exceptions to the foregoing assumption in that the lifetimes of vibrational 
excited states at room temperature can be quite long. However, it is probable 
that such excitation is removed in a few collisions with the wall. Then the 
assumption of a negligible population of excited molecules is valid. Although 
the population of excited states may be small enough to have negligible effect on 
the solution of the diffusion problem, the excited molecules may exist long 
enough to permit ultraviolet absorption studies . 

Recombination is not the only chemical reaction possible with reactant 
species 1 and 2. In an experiment with oxygen, for example, recombination com- 
petes with ozone production. The chemical equation for this reaction is 


[0] + [0 2 ] + [M] [0 3 ] + [M] 


(58) 


Ozone formation Is, in turn, followed by the reaction 


[O3] + [0] -► 2[0 2 ] 


(59) 


A complete "source" or "sink" term will consist of a linear sum of terms each 
representing the rate of loss or gain, indicative of a possible chemical reaction, 
of a respective species . The problem then becomes quite complicated because of 
the number of unknown parameters that enter . By adjusting the amount of molecular 
oxygen available, that is, oH 2 , it is possible to control the predominance of 
either chemical reaction. It will be assumed that the parameters are such that 
recombination (eq. (57)) is the dominant reaction. The sink term associated with 
this reaction (ref. k, p. 8) is given by 


-Ki^nininj^ 


(60) 


for the rate of loss of species 1 due to a triple collision with reactants n x 
and a third body M. The constant K r y is called the reaction or recombination 
rate constant and has the dimension (cc/particle) 2 sec -1 . The subscript M 
indicates that the constant may depend on the kind of third body involved. 
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Let K ri , K ra , and Kr 3 be the rate constants associated with reactions 
resulting from a triple collision wherein the third body is species 1, 2 , and 3, 
respectively. Then the total rate constant for a change in species 1 resulting 
from collisions is 


Ki = -(K ri nx 3 + Kp^n^ng + Kj» 3 nx 2 n 3 ) (6l) 

By expansion in terms of A using the pertinent relations in equations (ll), 
the following relations are obtained 


oKi = 0 (62a) 

iKx = -[ (K r£ A) inj 2 o n 2 + (^r 3 A) x 11 ! 2 0 ^ 3 ] (62b) 


In equation (62b) , the quantities (Kr 2 A) and (Kr 3 A) are taken to be of order 
zero with respect to the expansion for small A. This is done in order to make 
the sink term of the same order as the diffusion terms. In the experiment the 
parameters will be chosen in just this way; that is, the pressure and the geom- 
etry of the diffusion tube will be adjusted so that the homogeneous recombination 
and the diffusion terms are of comparable magnitude. 

When equation (o2b) is substituted into equation (25e) and the number 
density expressions are made dimensionless (using eq. (25&))j bhe resulting 
expression becomes 


x k x = - 


ch 2 AK r2 ^ oM- 2 + ( ~ o n2 ^r 3 ) ihi 2 oM-3 


D 


(63) 


Letting k r2 and k r3 be the dimensionless constants associated with the reaction 
constants Kr 2 and Kp 3 , where 


k r2 o n2 A J K r2 


(64a) 


kr 3 “ ( 3" on 2 A j K r3 


(64b) 


and 


i^x 


= -(kr 2 od 2 + ^r 3 0^3) idi 2 


(65) 
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the zero-order term with respect to our parameter fl using equation (31b) 
becomes 

1^1 = “(kr 2 0^2 + ^r 3 0^3) (1^1) (66) 

If k r2 = kp 3 , equation (66) then becomes 

1^1 = -k r2 (x^i) (66a) 


For the experiment discussed in the introduction, solutions based on a sink 
term defined by equations (6l) through (66a) are needed. This requires solution 
of a nonlinear equation. A numerical procedure for solving this equation by 
machine computation is being developed. If only a qualitative understanding is 
desired, the equations can be linearized since closed-form solutions can be 
found for the lin earized probl em. Equation (66a) can be linearized by replacing 
(l^i) 2 with iPi(il^i)^ where ( xPi ) is a constant average value. This method is 
used, rather than a conventional linearization about a constant value of 1P1, in 
order to obtain a zero value of Jkx at = 0, in qualitative agreement with 

the nonlinear expression. Then equation (66a) is replaced by 


xkx - -k r xM-i (67a) 

where 

k r = ~kr 2 SI = " \ kp 2 (67b) 


Here Jpi is taken equal to l /2 because in the results to be given £px varies 
between 0 and 1. 

Equation (67a) can also be interpreted as a result corresponding to the 
reactions represented by equations (58) and ( 59 )> as follows. If it is assumed 
that any ozone molecules formed are quickly converted to oxygen molecules by the 
reaction described by equation (59) j then the total sink term associated with 
equations (58) and ( 59 ) is linear in [ 0 ], That is, it is assumed that in a very 
short time (compared to diffusion time) after an ozone molecule is formed by the 
chemical reaction (eq. (58)), the molecule will undergo a binary collision and 
that the efficiency of reaction is so high that the collision will cause a 
chemical change. This reaction is thus presumed to depend only on the species 
number density nx . With these assumptions our general system can still be 
treated as a three -component system, but the sink term can be treated, with 
respect to species 1, as linearly proportional to the concentration of species 1. 
In terms of these assumptions the sink term with respect to species 1 corre- 
sponding to ozone formation (eq. (59)) is 
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( 69 ) 


The term corresponding to recombination (eq. (58)) is 


-K r nin 2 o n 


Here differences that might arise because of the kind of third body involved in 
the triple collision required for equation (58) are not considered • When equa- 
tions (68) and (69) are summed, the expression obtained represents the loss of 
species 1 in the assumed chain of reactions; that is, the sink term becomes 


Kjl = -(Kp ? n^ + K r nin 2 


( 70 ) 


When a procedure similar to that leading to equation ( 70) used, equation (63) 
then becomes 


Letting 



iMa + 


L i K r on 2 


lM-l oM-h 


k r 


l 2, 

n t r J 


D 


K r 



( 71 ) 


( 72 ) 


Equation (71) can be represented by the expression 


ik x = -k r ini 


( 73 ) 


In s umma rizing we observe that a quantitative representation of the sink 
term is given by equation (66) or (66a) , but the partial differential equations 
(eqs . ( 54 a), (54b), and ( 54 c)) resulting from this substitution require electronic 
machine computation because of their nonlinearity. However, qualitative results 
may be obtained from the closed-form solutions resulting from the linear 
approximation to the sink term, equation (67a) or (73) - 


Boundary Conditions (including Wall Reaction Effects) 


The equations have been kept general up to this point since their derivation 
is not affected by any specific choice of coordinate system. In the discussion 
of the boundary conditions a specific configuration will be used, as described 
in the introduction. The experimental technique motivating this analysis utilizes 



a flowing gas made up of known amounts of an inert gas and molecular oxygen. 

The oxygen is partially dissociated in such a way that there is a periodic vari- 
ation in the concentrations of both the atoms and the molecules. In this respect 
the concentration of the various species at the base of a cylindrical diffusion 
tube, where the recombination process is being studied and which Is tangent to 
the flowing gas, is expressed in terms of the total gas density by 


n(o) = ni(o) + n 2 (o) + n 3 (o) 


( 7 ^) 


The zero in parenthesis signifies that these concentrations are evaluated at the 
inlet end of the diffusion tube . 

When equation (74) is expanded about a uniform state in terms of the 
parameter A, as was done to obtain equations (lla) through (llg) and when 
coefficients of like powers of A are equated, the following equations are 
obtained 


0 n( 0 ) = Q n 2 (o) 

+ o n 3 ( 0 ) 

(75a) 

in(o) = ini(o) 

+ 1 ^ 12 ( 0 ) + 1 ^ 3 ( 0 ) 

(75b) 


The zero-order terms are constant . Also, as was described earlier, there are no 
external variations in pressure. The following identities result: 


Q n(o) = Q n (76a) 

on 2 ( o ) = o n 2 ( 7bb ) 

on 3 (o) = (76c) 

ini(o) + in 2 (o) + fn 3 (o) = 0 (76d) 


To determine the boundary conditions on the inlet surface of the diffusion 
tube, expressions must be found for ini(o), in 2 (o), and jn 3 (o) . Before these 
expressions can be determined, it Is necessary to introduce a dimensionless 
function that represents the time variation applied to the dissociation process. 
Let fi(t) be a periodic function representing the dissociation time function for 
species 1, such that 


ni(o) = n 1 (o) max f 1 (t) 


(77) 
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represents the concentration of species 1 on a two-dimensional surface repre- 
senting the inlet boundary of the diffusion tube. The dimensions of the system 
are such that ni(o) max can be considered constant with respect to the radial 
variable on the base of the diffusion tube. If the dissociation is periodic, the 
dissociation time function can be expressed as 

f l( t ) = 1 - COS Wt ( 78 ) 

where w is the frequency of the time variation. 

When equation ( 77 ) is expanded in terms of the parameters A and ft and the 
density is made dimensionless, it then follows that 

2m(o) = fi(t) (79) 


Equation ( 79 ) is the inlet boundary condition applied to equation (49a). 

The variable in 2 (o) is related to ini(o) through the dissociation process. 
This is expressed by the equation 


mi ini(o) + life in 2 (o) = 0 


( 80 ) 


This relation neglects diffusion processes in the outer flow tube . Making 
equation (80) dimensionless and then substituting equation (79) into the resulting 
relation we have 


iha(o) 



( 8 la) 


When the relations for ihi(o) and °m- 2 (o) and equation (76d) are used, the 
equivalent condition for species 3 is obtained. It then follows that 



(8lb) 


Sufficient information is now available to determine the boundary conditions 
on the inlet end of the diffusion tube for use in solving equation (49b). The 
condition for equation (49b) is determined by substituting equations ( 8 la) and 
( 8 lb) into equation (49c). The result is 


v(o) = - | fi(t)(l + g) 


( 82 ) 
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It is now necessary to determine the effect of the other surfaces that 
constitute the boundaries of the side-arm diffusion tube. Here the normal gradi- 
ent in the number density of all species is related to a fraction 7 , which 
represents the fraction of atoms that recombine upon collision with a surface. 

For this problem only the case where the stream velocity is zero will be consid- 
ered. The fraction 7 is defined in the following manner. If J in represents 
the effusive flux of particles of species 1 into a surface, and J£ n represents 
the effusive flux in the opposite direction, the net flux of atoms in the region 
of the surface Is given by 

(■In - Jin = /Jin) at vail ( 8 3 ) 

The effusive flux which represents the number of particles that cross the 

plane ds per unit time per unit area in the positive direction with respect to 
ds, is determined from (see ref. P- ^56) 


T + - ds 

~ — 

ds 



fiVx dV x 


(84) 


The region of integration S includes only those values of V in which have a 
positive component with respect to the normal to the plane ds. The distribution 
function is given by (ref. 3, ch. 7) 




1 + n • djCiJCWi 2 ) 

j 


(85) 


where fx is the Maxwellian distribution given by 

,3/2 


r o _ _ I m i 
1 1 ' 2rtkT 


•mxV x 2 / 2 kT 


(85a) 


and CxJ is a^scalar parameter which is determined in reference 3, chapter J. 
The quantity dj is given by 


d = A (Hi + (li - n H ^ Zn p 

*) Sr \ n / \ n P ) dr 


and the reduced velocity Wj_ 


is defined by 


Wi = 


pi Vi 

v/2kT 1 


( 86 ) 


(87) 
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When equation (85) is substituted into equation ( 84 ), in suffix notation 
(a = 1, 2, 3) with the usual convention that repeated dummy subscripts denote 
summation, it follows that 


JiZ = 


fi°VxZ 


dV] 


“I 




( 88 ) 


where the integrals are defined by 

ifja = P [ WiaWxZfx^i^Wi 2 ) dV x (89) 

The first integral on the right side of equation (88) when integrated becomes 


r n i v i 

4 


(90) 


where 


vi = 



(90a) 


is defined as the mean molecular speed. From symmetry properties alone and the 
fact that the integrations are over the entire space of Vix and V ly the only 

nonvanishing integral in equation (89) is ifjZ- The integrand of this integral 
is even in the reduced velocity W x and is therefore equivalent to 


1 

6 


W x ^f x ^C X ') dV x 


(91) 


Noting that this integral is the same as that occurring in the definition for the 
multicomponent diffusion coefficients (ref. 3> P- ^ 79 ) > we obtain the identity 


lijz - \ =1 Eij (92) 

When equations (90) and (92) are substituted into equation (88), the effusive 
flux becomes 

jiz = ^ n i^i + (93) 

j 
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By using a similar method of analysis we obtain for the number density flux in 
the opposite relative direction 


JiZ = p nivj. 
4 



djzDij 


( 94 ) 


Equations ( 93 ) and ( 9*0 are now substituted into equation (83) and after making 
the resulting expression dimensionless and expanding it in terms of A and it 
follows that 


o £ 
oP 



3 


ojDxj 8 

D sT 



evaluated 
at wall 


(95) 


Upon substitution of the multicomponent diffusion constant relations, equa- 
tions (21a) and (21b), equation (95) reduces to 


5 



(2m) 


ds 

ds 


■iM-i^Wx 


J at wall 


where 


(96) 


&Wi 


DD23 f 11A 7 

aD 12 B 13 \ 4 D ) 7 

X ~ 2 


(96a) 


Equations (96) and (79) are sufficient to define the boundary conditions for 
equation ( 49 a). It only remains now to determine a relation similar to equa- 
tion (96) as a boundary condition for the quantity v appearing in equation ( 49 b). 
Since species 3 is chemically inert, it does not react upon collision with the 
wall. The surface equation equivalent to equation (83) for species 3 is 


(<J 3 n _ J 3 n)at wall - 0 


( 97 ) 


The number fluxes J3X1 and J3n are determined by the same scheme used to derive 
equations (93) and (9*0 t o obtain 



(98) 
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When operations similar to those utilized in going from equation (95) to (96) 
are used, it follows that 


f <**«> • £ - ^ 


(l - Oh2c) 


ms 


oM -3 ~ (l - d) + 0^2(1 - c) 




(99) 


at wall 


after equations (96) and (99) are substituted into an expression found by taking 
the gradient of equation (76), the corresponding equation for species 3 is 
obtained. The result is 


• 








4 (&* 

1 

0 

~iMa • 

OM-3 

a (1 . d) - oto _ 

kwi 


(100) 

ds 


(l - oM- 2 c ) 


V 





✓ 

at wall 


for v 

becomes 

fbv 

\W ' 

ds , 
ds 

= lhlkw 2 ) 

/ at wall 



(101) 


where 


^w 2 



mi 

m 2 


+ 




( 101a) 


When the boundary equations (96) and (lOl) are applied to a three- 
dimensional cylindrical diffusion tube, there will be a separate equation for 
each surface (cylinder and end wall). In the event that the various surfaces 
are constructed of different materials, the fraction v will be different and 
will depend on the material used • ( See ref . 5 • ) 


DISCUSSION 


The proposed experiment which prompted the analysis contained in this 
report will utilize the decay of the time variation in species concentration in 
determining the diffusion and homogeneous recombination coefficients of oxygen 
atoms. The preceding analysis establishes the mathematical formulation for 
this problem. A complete solution to the equations that constitute a formulation 
of the problem, that is, the determination of the concentrations of all species 
for all coordinates in space and time, will not be covered in this report. 
Instead, a solution will be given that illustrates the dependence of the concen- 
tration variation decay on chemical recombination and diffusion. This solution 
will apply to a one -dimensional system where diffusion and recombination occur 
with zero stream velocity in a region between two infinite planes and where the 
species concentration is being varied on one plane. At the other plane the 


28 



condition of equation ( 96 ) is imposed. In this example the linear sink term is 
used and heterogeneous recombination is neglected. The time variation is that 
given by equation ( 78 ) • 

The solution is found by following a scheme similar to that given by 
Sommerfeld (ref. 6) in his solution to the diurnal temperature problem. It 
follows that a solution which includes parameters that are determined so that 
the differential equation and the boundary conditions, equations ( 78 ) and ( 98 ), 
are satisfied, is given by 


= u 0 (|) + Real [exp(i¥T)u x (|) ] 


( 102 ) 


where 


u^o) =1 Z = 0,1 


(l02a) 


When equation (102) is substituted into equation (54a), where the sink term is 
given by equation ( 67 a) or ( 73 ) j the resulting expression is a linear ordinary 
differential equation in U2- The solution of this equation is 

u 1 = A p+ exp( P + £ ) + Ap_ exp(P_£ ) (IO 3 ) 

where (P + ) and (P_) are roots of the characteristic equation 

P 2 - (l 4- iWZ) =0 Z = 0,1 (104) 


and where Ap + and Ap_ are determined by the boundary conditions, equations ( 96 ) 
and (l02a). When the wall reaction is neglected, the solution becomes 


o. \ cosh L/k^ (| - l) _ n 

iPiU,t) = - ^ ^ + Real 

cosh Jkp 


a iWT cosh [p x (£ - l)] 
cosh Pi 


( 105 ) 



This solution consists of one time -independent term superimposed on a term that 
includes time dependence. The average concentration is given by the time- 
independent term, that is, the first term on the right side of equation ( 105 ) . 
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The concentration time variation is given by the second term. To illustrate the 
dependence of the decay on diffusion and recombination, a depth parameter £ x 
will be introduced that represents a distance at which the variation in concen- 
tration is reduced to e - - 1 - of its initial value . The depth parameter is defined 

by 


cosh [Pi(|i - l). 
cosh Pi 


e 


-i 


(106) 


where the expression on the left of equation (106) is the magnitude of the time 
variation that results from the second term on the right side of equation (105) . 

Figure 2 illustrates the dependence of the depth parameter | 1 on frequency 
nnd sink strength. It is to be observed that for high frequencies the decay in 
concentration occurs close to the inlet plane and is little affected by the sink 
strength k r . This follows if it is noted that when the dimensionless frequency 
W/2rt is of order 20 or greater, the l/e decay characterized by li occurs 
within l/lO of the diffusion tube length and depends very little on the value 
of the sink term parameter k r . As the frequency is reduced, there is a greater 
dependence of the depth parameter on the sink strength, and the decay occurs 
farther from the inlet . Also, in the limit of large sink strength the decay 
occurs near the inlet plane and is little affected by the frequency. For 
instance, we note that when k r is of order 100, a change in dimensionless fre- 
quency from 0 to 30 produces a change in the decay depth from 0.10 to 0.08.^ 

From the preceding discussion we observe that to obtain maximum resolution in 
the determination of k r , it is desirable that the frequency of variation be low 
but not zero since the limiting case of zero frequency (or alternatively, the 
steady state problem) is not of interest here. The reason for this is that the 
variation in species concentration is itself a technique for making measurements 
on the reaction rates. The appropriate frequency will, in part, depend on the 
electric filter circuitry efficiency in separating the time varying signal from 
the background. The detection of the concentration variation may also be 
improved, particularly for large kr, by superimposing the variation in concen- 
tration on a steady gas stream of known velocity. Such a scheme serves to trans- 
port the atoms from their source to the region of observation in a shorter time. 
This is necessary in the event that the recombination time is small and the 
recombination occurs before the atoms are in the region where the .observations 
are being made (this case is characterized by large k r ) . The utility of the 
proposed experiment arises out of the sensitivity available in photometric 
detection techniques associated with spectral absorption experiments for meas- 
uring relative changes in concentration. It is therefore expected that the 
decay in the concentration variation may be measured with great accuracy. 


CONCLUDING REMARKS 


A set of nonlinear and coupled fundamental equations describing diffusion 
processes in general, obtained from rigorous kinetic theory, has been used, to 
derive a set of simplified expressions applicable to an unsteady problem in 
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which the species concentrations are made to vary in a prescribed periodic 
manner. With the set of simplified expressions the concentrations of each 
species in a ternary gas mixture can be determined as functions of time and 
spatial coordinate. It was necessary to complement the derivations by introduc- 
ing chemical kinetics in discussing one term in the fundamental equations (the 
term representing the recombination processes ), and by introducing wall reaction 
effects to establish the boundary conditions necessary for the final set of 
simplified expressions. The term representing the recombination processes intro- 
duced a nonlinearity in the simplified equations. However, a simplifying assump- 
tion could be made which allowed for closed-form solutions to facilitate gaining 
a qualitative understanding of the unsteady problem. A solution based on this 
approximation was given . It required the same methods that were applied in 
determining the diurnal or annual temperature variation of the earth. In fact, 
there is an obvious analogy between the experimental problem mentioned above and 
the earth ? s temperature problem. In the earth T s temperature problem the vari- 
ation in temperature due to the daily solar fluctuation is confined to a thin 
surface layer. The experiment utilizes measurements on the "width 11 of the 
analogous concentration layer in determining the recombination rates. The 
qualitative picture gained from the solution illustrates that care must be taken 
in choosing the parameters of the problem. For example, the technique is appli- 
cable only if the variation frequency is sufficiently low that the effect of 
chemical reaction can be observed; also, the experimental method may require 
that the concentration be superimposed on a gas stream should the reaction rates 
be too large and the decay confined too close to the inlet boundary. 

The next logical development, measuring the recombination and diffusion 
rates, requires the determination of the solution of the equations with the 
nonlinear recombination term. This solution will, of necessity, involve recourse 
to numerical methods . Such a program is under way, and preliminary exploratory 
work indicates that there will be no difficulty in achieving a solution to this 
more difficult problem. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Nov. 9, 1962 
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Figure 2.- Curves relating the depth parameter, a parameter which characterizes 
l/e decay in concentration variation, to the dimensionless frequency for 
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various values of the linear sink term coefficient k r - 
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